Wigner crystallization in the two electron quantum dot 

A. Matulis 

Institute of Semiconductor Physics, Gostauto 11, 2600 Vilnius, Lithuania 

F. M. Peetersl 

Departement Natuurkunde, Universiteit Antwerpen (UIA), Universiteitsplein 1, B-2610 Antwerp, Belgium 

(February 1, 2008) 



Wigner crystallization can be induced in a quantum dot by 
increasing the effective electron-electron interaction through 
a decrease of the electron density or by the application of a 
strong magnetic field. We show that the ground state in both 
cases is very similar but the energy scales are very different 
and therefore also the dynamics. 
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I. INTRODUCTION 

During the last ten years quantum dots have attracted 
a lot of interest both experimentally and theoretically [Q . 
Recently, much attention was paid to the investigation 
of quantum dots in strong perpendicular magnetic fields 
where a rich structure of cusps and steps in the chemi- 
cal potential n(N, B) (as a function of confined electron 
number N and applied magnetic field strength B) was 
observed The above structure is related to changes 
in the ground state electron density in the presence of 
a magnetic field, and it initiated numerous theoretical 
investigations after Chamon and Wen || proposed the 
quantum dot edge reconstruction. The electron density 
in a quantum dot is the result of the interplay of the 
repulsive character of the electron interaction and the 
attractive forces caused by the confining potential, the 
magnetic field and the electron exchange interaction. If 
the magnetic field is strong enough, the overlap between 
the electron wave functions becomes less, the electron in- 
teraction will dominate, and a ring of electrons at the dot 
edge is formed. It was also found that if the magnetic 
field is increased further the above electron ring becomes 
unstable, and a ground state with a broken rotational 
symmetry appears. The possibility of the appearance of 
spin waves Jl[ and charge density waves |5|,|| have also 
been reported. 

A quantum dot with a density profile consisting of rings 
with lumps reminds us to the Wigner crystal which 
is the ground state of the classical electron system in a dot 
Pj in the absence of a magnetic field. In the latter case 
the Wigner crystal occurs when the potential energy (the 
inter-electron interaction and the confinement potential) 
dominates over the electron kinetic energy. This is just 
the classical limit in a quantum problem. Therefore, clas- 
sical or quasi-classical methods should be adequate for 



the description of the Wigner crystal. But from a first 
sight it appears that such a classical limit is not reached 
when a strong magnetic field is applied. The above men- 
tioned electron density reconstruction was revealed when 
the magnetic field is so strong that the electrons occupy 
the lowest Landau level. Thus, the electron motion quan- 
tization is essential, and consequently, the kinetic energy 
exceeds the potential energy due to the Coulomb interac- 
tion. Nevertheless, due to the degeneracy of the Landau 
level this large electron kinetic energy is actually frozen 
out, and Wigner crystallization results from the same po- 
tential energy as in the classical case without magnetic 
field. Actually the above crystallization is a result of 
different energy scales in the electron system under con- 
sideration. 

The purpose of the present paper is, by using an ex- 
actly solvable model of two electrons in a dot, to illus- 
trate the conditions under which a Wigner crystal can be 
formed in the case of a strong magnetic field (B), and to 
compare this quantum Wigner crystal with the classical 
zero magnetic field one. The paper is organized as fol- 
lows. In Sect. II we present our model. Sect. Ill gives 
an introduction to the quasi-classical approach for the 
B = case. Sect. IV discusses the general B ^ case. 
Our conclusions are presented in Sect. V. 



II. MODEL 

We consider two electrons with effective masses m* 
moving in the z = plane and which are confined by 
a two-dimensional harmonic potential of characteristic 
frequency loq. A magnetic field is applied in the z di- 
rection and described by the vector potential chosen in 
the symmetric gauge A = [B x r]/2. The corresponding 
Hamiltonian can be separated into two parts 



H = H, 



(i) 



(see, for instance |Uj,|llJ) which represents the center- 
of-mass and relative motion (with corresponding coordi- 
nates R = (ri+r2)/2 and r = ri — Y2). In dimensionless 
form those parts can be written as follows: 
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The system energy is measured in Tiloq units, and the 
coordinates in o,q — y / nJm*ujQ units. The symbol 
A = ao/a,B stands for the electron interaction coupling 
constant which is the ratio of the characteristic dot di- 
mension ao and the Bohr radius as — eh 2 /m*e 2 . The 
magnetic field strength B is measured in ^q/tto^ units 
where $o = ~xTic/e is the magnetic flux quantum. 

We will concentrate ourselves to the study of the coor- 
dinate wave function part of the ground singlet state 



*(ri,r 2 ) = *(R)V(r), 
the corresponding electron density 

p(r) = / d 2 r x f d 2 r 2 |*(n,r 2 )| 2 p(r) 



(3) 
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and the correlation function 



K(r,r') = J d 2 n J d 2 r 2 |*(n,r 2 ) 



x {p(r)p(r') - 6(r - r')p(r)} 
= 2|$({r + r'}/2)| 2 |^r-r')| 2 . (5) 

It is the latter function which enables one to determine 
whether, or not the system is in the Wigner crystal state. 

III. ZERO MAGNETIC FIELD CASE 

We consider first the simplest case in which no mag- 
netic field is applied. Then the center-of-mass motion 
part is trivial. It just represents the two dimensional 
harmonic oscillator motion which has no relation to the 
Wigner crystallization problem. Its ground state energy 
is Er = 1 with the corresponding wave function 
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Due to the cylinder symmetry the relative motion wave 
function part can be written as 



1 



V'(r) = exp(irmp) R(r) , 
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where the radial function R(r) is obtained by solving the 
one dimensional eigenvalue problem as determined by the 
Hamiltonian 
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As was already mentioned strong electron correlation 
(and the Wigner crystal as well) occurs when the poten- 
tial energy dominates over the electron kinetic energy, 
i. e. when A — > oo. In this interesting limit the eigenvalue 
problem is strongly simplified and can be solved by ana- 
lytical means. Indeed, in the case of A — ► oo the potential 
(pb| ) has a minimum close to the point ro = (2A) 1 / 3 (see 
solid curve in Fig. |). The potential can be expanded 
into a (r — ro) series (dashed curve in Fig. |l|) 



^A)^(2A) 2 /3 + ^-, ) 2 + ^ 



(9) 



which actually coincides with the quasi-classical expan- 
sion into a A~ 2 / 3 series ]Tjj| . The form of the above ex- 
pansion clearly indicates the energy scales of the different 
excitations in the quantum dot. The first term (A 2 / 3 ) 
is just the classical dot energy which can also be ob- 
tained using the hydrodynamic approximation |12[ ]. In 
the A — > oo limit this energy dominates. Thus, we are in 
the region of the classical Wigner crystal. The solution 
of the Schrodinger equation with the harmonic term in 
the expansion of the potential (|^) leads to A-independent 
equidistant electron ring vibration excitations (shown in 
Fig. by the horizontal dashed lines). Note that the 
separation of those vibration levels are of order 1, and 
consequently, much less than the classical potential en- 
ergy. The last term in the expansion (^J) describes the 
rotation energy. It leads to a small (of order ~ A~ 2 / 3 ) 
splitting of the rotation levels as is shown by the encircled 
part in Fig. [l] which is enlarged. The spectrum is thus 
similar to the one of molecules with bands of rotation 
levels, attached to each vibration level. 

For the expanded potential (^) we obtain the following 
quasi-classical radial ground state wave function, 



P(r) = 7Vexp{-a(r-r ) 2 }, 



(10) 



where a = \/3/4 and the normalization N — 

(2aA)VV /2 . 

Inserting the above relative motion wave function ( |l0| ) 
together with the center-of-mass motion function (^|) into 
expressions (D,||) we obtain the electron density 
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and the correlation function 
AN 2 



K(r,r') = ^-exp{-(r + r') 2 /2} 

x exp{-2a(|r-r'| -r ) 2 } . (12) 
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In expression (11), Io is the Bessel function and 7 = 
8a/(l + 4o) w 1.268. 

A schematic picture of the above functions is shown in 
Fig. H by the shadowed regions. In Fig. ||(a) the elec- 
tron density is depicted which is mainly concentrated 
on a thin ring. Notice that it has no lumps, and ex- 
hibits the cylindrical symmetry of the system Hamilto- 
nian. Wigner crystal state can be seen in the correlation 

function @0 which is P lotted in Fi g- @( b )- The latter 
is actually the same as the conditional probability distri- 
bution. One electron is fixed in position r' (indicated by a 
cross in Fig. ||(b)) and the density of the other electron is 
then mainly concentrated in the opposite position which 
clearly indicates the presence of Wigner crystallization in 
this quasi-classical limit. 



IV. STRONG MAGNETIC FIELD CASE 

Due to the symmetric gauge which preserves the cylin- 
drical symmetry, both (i. e. center-of-mass and relative) 
wave functions can be written as a product of the or- 
bital exponent and the radial part as in expression (0). 
The corresponding radial Hamiltonians (for center-of- 
mass and relative motions) can be presented as follows: 



has energy Er — E r = 1/2 labeled by integer posi- 
tive momentum (M and to) values. Mathematically this 
degeneracy is a consequence of the equivalence of the 
zero center-of-mass Hamiltonian potential (M/R— R) 2 /A 
which is shown in Fig. ^|(a) where the orbital momentum 
is indicated by the corresponding integers. For any mo- 
mentum M the potential curve has a minimum equal to 
zero at the position iZ m i n = ^/M. A similar potential is 
obtained for the zero order relative motion equation, and 
will therefore not be discussed. 

Next we take the perturbation terms into account. For 
the case of center-of-mass motion the potential is shown 
in Fig. |^(b). The potential is composed of the same 
curves as in the case without perturbation, but they are 
now moved up by the amount R 2 /B 2 shown in Fig. ||(b) 
by the dotted curve. The lowest minimum is obtained 
for M = which implies that the center-of-mass motion 
wave function is the same as in the case without magnetic 
field (§). 

This is not the case for the relative motion where all 
potential curves are shifted up by the amount V r (r, B) 
as s hown in Fig. ||(c) by the dotted curve. According to 
Eq. (16b) this curve increases for both r — > and r — > 00 
and reaches a minimum value at 
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Because we are interested in the asymptotic limit B — ► 
00, it is convenient to scale the variables and the Hamil- 
tonian as follows 



vB 



-1/2 



R->RB- 1/2 , H^HB, 



(14) 



in order to have the expansion parameter B 1 explicitly 



in our problem. Inserting (14) into expression 
arrive at the following Hamiltonian: 
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where in the asymptotic region B 
Vr(R,B) = ^, 



the terms 
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can be treated as small perturbations. 

The solution of the Schrodinger equations for the zero 
order Hamiltonians (without the terms Vr and V r ) leads 
to the degenerate Landau levels where the ground state 



Minimizing the potentials with respect to the relative an- 
gular momentum for r — r m i n we obtain the correspond- 
ing orbital momentum 



2/3 



B. 



(18) 



Consequently, the relative motion wave function is 
peaked on the ring of radius r$. When the magnetic field 
strength increases the radius r tends to infinity, while 
the orbital number is growing as well. Expanding the 
potential with the lowest minimum in the vicinity of the 
equilibrium point ( p^7|) we find 
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which means that the thickness of the ring remains con- 
stant. 

The layout of the energy spectrum in the large mag- 
netic field limit is shown in Fig. |^. Notice that the spec- 
trum has two different energy scales as in the B = 
case. But the physical meaning of those scales is quite 
different. The largest energy scale is the electron kinetic 
energy. And the rotation levels are split due to the in- 
terplay of the Coulomb interaction between the electrons 
and the confinement potential. 

Going back to our original units before the scaling jl^ ) 
we find that the ring radius is r m ; n = (2A) 1 / 3 , which is 
identical to the zero magnetic field case. But now the 
thickness of the ring is ~ B^ 1 / 2 , and it tends to zero as 
the magnetic field strength approaches infinity. 
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Surprisingly, we arrived at the same situation as was 
found for the case without magnetic field. The relative 
motion wave function is located on a ring whose diameter 
greatly exceeds its thickness. Thus the conclusions of 
previous section concerning Wigner crystallization in the 
quantum dot are also valid in the strong magnetic field 
case. While for B = the Wigner crystal state is reached 
for A 3> 1 we find that a sufficiently strong magnetic field 
can crystallize the system for any A 7^ 0. 



V. CONCLUSIONS 

The occurrence of Wigner crystallization depends on 
the ratio of the distance between the electrons (7) and 
the characteristic dimension of the single electron wave 
packet (a) . Actually this ratio is a measure of the electron 
density in the quantum dot. In our case I is given by the 
radius r min of the ring in the correlation function, and a 
coincides with its thickness. This ratio is 
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Here Ib = \Jfic/eB is the magnetic length. 

The larger this ratio the more pronounced the Wigner 
crystal state is. Thus, both a strong magnetic field and 
strong interaction favors Wigner crystallization, but in a 
different way. The electron interaction makes the system 
less dense by enlarging the inter-particle distance. The 
magnetic field also makes the system effectively less dense 
but by compressing the single electron wave packages. 

Also we would like to note that mathematically in both 
cases the static Wigner crystal properties can be calcu- 
lated by the same method, namely by the minimization 
of the classical potential energy which is composed of 
the electron interaction and the confinement potential. 
Nevertheless the physical meaning of that calculation is 
quite different. In the strong electron interaction case the 
potential energy dominates making the whole problem 
quasi-classical, while in the case of strong magnetic fields, 
the problem is essentially quantum mechanical (the Lan- 
dau level energy is dominating). But due to the degen- 
eracy of the problem the system is guided by the same 
potential energy as for the B = case but with different 
energy scales. Therefore, although the Wigner crystal 
configuration is the same, one can expect different dy- 
namics in the two limiting cases. 
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FIG. 1. The potential energy for the relative motion (solid 
curve) and its quadratic expansion (dashed curve) around its 
minimum (j^) . The horizontal lines are the different vibration 
levels in units of huio- 

FIG. 2. Schematic view of the electron density (a) and the 
electron correlation function (b) in the quasi-classical limit. 

FIG. 3. The potential energy curves for different angular 
momentum in case of: (a) the non-interacting problem, (b) 
the center-of-mass motion including confinement, and (c) the 
relative motion including confinement and electron-electron 
interaction. 

FIG. 4. Energy spectrum in the limit of a strong magnetic 
field. The energy is given in units of hui c . 
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